4532 


Predicting Contaminant Removai 
Effectiveness of Three Air Distribution 
Systems by CFD Modeling 

Daniel J. Hirnikel, P.E. Peter J. Lipowicz, Ph.D. Raymond W. Lau, Ph.D, 

Member ASHRAE Member ASHRAE Associate Member ASHRAE 


ABSTRACT 

This study compares the. contaminant removal effective¬ 
ness (CRE) of three air distribution systems for a bar/restau¬ 
rant setting by use of computational fluid dynamics (CFD) 
modeling. The air supply distribution and exhaust arrange¬ 
ment were modeled for a mixed air system, a directional air 
flow system (where air moves in a unidirectional flow pattern 
across the space,from the high-pressure supply area to the low- 
pressure exhaust area), and a displacement ventilation system 
{where supply air is delivered at low velocity close to the floor, 
allowing thermal plumes to be created). The CRE of each space 
was determined for both particulate and carbon monoxide 
dispersions under two different ventilation rates. A commer¬ 
cial CFD software package wai used to model a hypothetical 
restaurant with bar in a steady-state condition with the inclu¬ 
sion of objects, such as a bar and multiple tables, heat sources 
such as people standing and seated, pollutant sources, multiple 
four-way supply air diffusers, and an exhaust air grille. The 
hospitality spaces’ smoking and nonsmoking zones had no 
physical separation. The results demonstrated that directional 
airflow systems reduce nonsmokers’ exposure to contaminants 
better than mixed air systems, and displacement ventilation 
systems reduce nonsmokers’ exposure better than both direc¬ 
tional airflow and mixed air systems. 

INTRODUCTION 

The purpose of any ventilation system is to provide a 
sufficient amount of outdoor air to remove contaminants and 
to maintain the interior environment at a comfortable temper¬ 
ature and relative humidity. The effectiveness of a ventilation 
system to perform these functions can be quantified in a 
number of ways e.g., air change efficiency of a system looks 


at air flow patterns only (Chow and Fung 1996; Nielsen 1993; 
Roos 1999), whereas contaminant removal effectiveness 
(CRE) quantifies tlie dilution andiemoval of acontaminant by 
a system. CRE is defined as the ratio between the steady-state 
concentration of the contaminant at the exhaust and the room 
or zone average value (Liddament 1993; Yagiou and With- 
eridge 1937)- 

The evaluation of CRE of three air systems was carried 
out using computational fluid dynamics (CFD) modeling. 
CFD modeling software has only recently become commer¬ 
cially available to heating, ventilating, and air-conditioning 
(HVAC) engineers in a package simple enough to learn and 
apply in a fairly short time. One virtue of CFD modeling is that 
it allows specific entry of details of the room that have rele¬ 
vance to the airflow. This feature is one that sets it apart from 
simple mass balance type models (Klepeis et al. 1996; Ott 
1999) in whicharoomislargelydescribedbyits volume. CFD 
models have been used recently to study indoor air quality 
(lAQ) problems, pollutant dispersions, and performance of 
HVAC systems (Chow and Fung 1996; Emmerich 1997; 
Gadgii et al. 1999), Indoor environments studied using CFD 
modeling included offices (Reynolds and Hedge 1999), 
research laboratories (Memarzadeh 1999), an operating 
theater (Tinker and Roberts 1999), and an industrial workshop 
(Peng and Davidson 1998). Previous studies involving the 
evaluation or comparison of different air distribution systems 
have been published (Jiang et al. 1992; Holmberg et al. 1999; 
Kolokotronj et al. 1995; Yeoh and Li 1998). Ventilation plays 
an important role for reducing environmental tobacco smoke 
(ETS) in tlie design for smoking areas (ASHRAE 1981,1989; 
Bohanon et al. 1998). The effect of ventilation on reducing 
ETS has been quantified in both laboratory studies (Cain et al. 
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(a) Mixed Air 



prevent smoke from entering nonsmoking areas or enhance 
the upward movement of smoke. 

The objective of this study was to compare the perfor¬ 
mance of three air distribution systems according to their abil¬ 
ity to remove contaminants, as estimated by their CRE 
(Enbom et al. 2000; Koganei et al. 1993; Liddament 1993; 
Roos 1999; Yaglou and Witheridge 1937). We quantified tire 
effect of ventilation under well-defined conditions of source 
strength, supply and exhaust size and location, and space 
geometry. The CFD model used is an ideal tool because it 
allows for the variation of ventilation design as an independent 
variable across different alternatives. Such variation is expen¬ 
sive in the laboratory and infeasible for a field study. 

In this study, a hypothetical hospitality space (arestaurant 
witli a bar) was modeled with numerous airflow obstructions, 
heat sources, and tobacco smoke-generating sources placed 
widiin the smoking section. Two components of ETS were 
modeled; particulate and carbon monoxide (CO). The restau¬ 
rant/bar was ventilated at two different airflow rates using 
three different air distribution systems, namely, a mixed air, a 
directional airflow, and a displacement ventilation system. 
The resulting distributions of the particulate and CO concen¬ 
trations were then obtained and used for the derivation of the 
respective CRE values. 

METHODS 

Computational Fluid Dynamics Program 


(b) Directional Airflow 



Figure 1 Room geometry in plan view. 

1983; Cain et al. 1987; Enbom et al. 2000; Ounnarsen and 
Fanger 1991; Kriihne and Pitzner 2000; Straub et al. 1993; 
Walker et al. 1997; Yaglou 1955) and some field studies 
(Bayer and Fischer 1997; Hedge et al. 1990; Hyvarinen et al. 
1997). Most design guides principally rely on the ability of 
ventilation air to dilute ETS. In practice, ventilation may have 
a lesser effect, due to poor supply air distribution, or a greater 
one, if it is also used to establish pressure gradients that 
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This study used a commercially available CFD program 
(Flomerics 1999) that could be readily applied by HVAC 
design engineers. This particular CFD program allowed a 
three-dimensional simulation to be carried out in a steady- 
state condition. The modeled space was divided into approx¬ 
imately 100,000 volume cells. The CFD program estimated 
variables such as pressure, velocity, temperature, and air 
constituent concentration for each cell, throughout the entire 
space in accordance with mass and concentration conserva¬ 
tion equations. 

The Model 

The room modeled was 10 m long, 7 m wide, and 3 m high 
(Figwre 1), consisting of a bar area (smoking area) and a 
restaurant area (nonsmoking area). The bar had nine people 
and one bartender standing around it. Each person generated 
approximately 100 watts of heat (ASHRAE 1997). 

One tobacco smoke source was placed on the bar 0.3 m 
away from the nonsmoking area. This location was 
purposely selected to be close to the nonsmoking section to 
evaluate the potential movement of particulate and the gas 
phase component, CO, into the nonsmoking region. The 
other smoke source was placed on a tabletop closest to the 
non-smoking area. The use of two tobacco smoke sources in 
the model was based on tire consideration of occupant load¬ 
ing of the space. It was assumed that each cigarette would 
last approximately eight minutes. Therefore, two cigarettes 
burning continuously were equivalent to fourteen cigarettes 
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per hour. A particulate emission rate of 1 rag/min for each 
cigarette was chosen. This falls into a range of emission rates 
used by other studies (Klepeis et al. 1996; Miller-Leiden and 
Nazaroff 1996). Two 1 mg/min particulate source would 
therefore generate 120 mg/h. Likewise, CO generation was 
modeled as 9 mg/min, which falls in a range identified by 
others (Klepeis et al. 1996; Jenkins et al. 2000). 

In accordance with the ASHRAE ventilation rate proce¬ 
dure (ASHRAE 1999), a bar with 18 people would require 
255.6 L/s of ventilation air, while the dining area with 23 
people wouldrequlre2i6.2L/s. As aresiUt, the modeled space 
required a total outdoor ventilation rate of 472 L/s. 

To investigate the effect of reduced ventilation on the 
CRE of the three air distribution systems, simulations were 
also carried out with 50% of the required ASHRAE ventilation 
rate, i.e., 236 L/s. 

Each tobacco source was modeled as 6 watts of heat 
(Waymack et al. 1997). This heat made the particulate and CO 
buoyant at the source. The diffusion coefficient used for CO 
was 2 X 10 ^ m^/s (Lundgren et al. 1979) and 2 x 10'^*^ m^/s 
(Lide et al. 1969) for the particulate, since these particles are 
0.2 micron in diameter (Jenkins et al. 2000). 

The exhaust grille was the location used to determine the 
whole space mean concentration for all cases modeled. Li a 
steady-state solution, all the mass flows into the space must 
equal all the mass flows exiting the space, which occurs at the 
exhaust grille. Therefore, the contaminant concentration in the 
exhaust air grille for the mixed air spaces, the directional 
airflow spaces, and the displacement ventilation systems was 
a constant (Tables 2 and 3) under the same ventilation rate. 

The spaces that were designed for mixed air and direc¬ 
tional airflow had ceiling-mounted foitr-way diffusers 
(Figures la, lb). The air was delivered into the room at WC 
with zero particulate and CO levels. This supply air was more 
dense than the room air; tlierefore, it dropped downward 
toward the floor due to gravity. There was one four-way 
diffuser located above the bar. It delivered 47.2 L/s, while each 
of three other ceiling diffusers delivered 141.6 L/s of ventila¬ 
tion air (a total of 472 L/s). Thus, each diffuser generated a 
primary air jet, which in turn set up multiple secondary air jets 
within the space. For the mixed air geometry, the supply air 
diffusers were located across the ceiling on both sides of tire 
space. For the directional air flow geometry, 90% of the supply 
air moved from the nonsmoking side of the space into the 
smoking side. 

In tire case of the displacement ventilation space, all the 
air was delivered to the space witli floor-mounted diffusers at 
low velocity (less than 0.1 m/s). All three cases modeled had 
their own unique airflow patterns, due to the heat sources and 
obstructions in each space. 

The CRE was determined by solving for the contaminant 
concentration in the exhaust grille divided by the mean 
concentration in the breathing zone of the nonsmoking region 
(Figures la, lb, and Ic). The face of a seated or standing 
person (between 1.2 m to 1.9 m elevation) in this nonsmoking 
section would be located within this specific region. 
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where CRE = contaminant removal effectiveness; Cj; = mean 
concentration in exhaust; Cg = mean concentration in supply 
air; Cg^ = mean concentration in occupied breathing zone. 

RESULTS AND DISCUSSION 

Mixed Air Systems 

The purpose of all mixed air systems is to create enough 
room air turbulence such that a uniform temperature and 
contaminant concentration exists. The CRE values (particu¬ 
late and CO) of the modeled restaurant/bar with the mixed air 
system at both ventilation rates were in a narrow range from 
1.0 to 1.2 (Table 1). The contaminant concentration was found 
to be typically within 15% of that at the exhaust grille (Table 
2). The higher ventilation rate resulted in a reduced 
nonsmoker exposure to contaminantsin mixed air systems due 
to dilution. Both particulate matter and CO are distributed 
similarly throughout the space, as is demonstrated by the CKE 
values. Both of these substances, although different in diffu¬ 
sion coefficient, mo ved as a result of the prevalent air currents 
within the space. 

The hospitality space was modeled with clean supply air. 
Tills was done to simplify tlie comparisons in the figures. 
Supply air is rarely found to be particulate free. However, as 
observed from formula 1, a contaminant in die supply air is 
subtracted out from both die numerator and denominator and 
results in the same CRE being determined as if it was clean 
supply air. 

Figures 2a and 3a display the concentration contour 
results for the mixed air systems at 472 L/s ventilation air. 
Figure 4a displays the particulate concentration contours for 
the lower ventilation rate. The concentration contours are 
shown as a plan cut through each space at 1.7 m elevation. The 
mixed air systems consistently (Table 1) exhibit the lowest 
CRE values in all cases. At the lower ventilation situation 
(Table 3) an uncomfortably high exhaust temperature (31“C) 
and higher contaminant concentrations were found. Also at 
the lower ventilation rate, less mixing took place widiin die 
space, resulting in CRE values slightly liigher than diose 
found at the higher ventilation rate. 

Directional Airflow 

The purpose of directional airflow systems is to supply 
clean air into the nonsmoking section and direct its movement 
into die smoking section in order to minimize the presence of 
smoke in the nonsmoking section. The CRE values (particu¬ 
late and CO) of die modeled restaurant/bar with the directional 
airflow system at both ventilation rates were consistently 
higher (Table 1} than the mixed air cases, indicating that 
contaminant concentrations in these nonsmoking regions 
were lower than those found when the mixed air systems were 
in use. 
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TABLE 1 

Contaminant Removal Effectiveness (CRE) Summary 


VentDatioii Rate 

Mixed Air 

Directional Airflow 

Displacement Ventilation 

472 L/s 

Particulate 

1.1 

2.7 

17.5 

CO 

1.0 

2.6 

15.8 

236 L/s 

Particulate 

1,2 

1,4 

7.4 

CO 

1.2 

1.4 

8.1 


TABLE 2 

472 L/s Ventilation Rate 


Parameters (unit of measure) 

Mixed Air 

Directional Airflow 

Displacement Ventilation 

Total supply airflow, L/s 

472 

472 

472 

Supply air particulate concentration, 

0 

0 

0 

Supply air temperature, "C 

14 

14 

16 

Exhaust flow, LVs 

472 

472 

472 

Exhaust temperature, “C 

23 

23 

25 

Particle concentration exhaust, 

70 

70 

70 

Pardcle concentration nonsmoking region, lig/m? 

63 

26 

4 

Particulate contaminant removal effectiveness (CRE) 

1.1 

2.7 

17.5 

CO concentration exhaust, pg/m^ (ppm) 

631 (0.6) 

629 (0.6) 

631 (0.6) 

CO concentration nonsmoking region, pg/m^ (ppm) 

624 (0.6) 

243 (0.2) 

40 (0.0) 

CO contaminant removal effectiveness (CRE) 

1.0 

2.6 

15.8 


Note: BoJd indicates cUlcuInted results from input data. 


TABLE 3 

236 L/s Ventilation Rate 


Parameters (unit of measure) 

Mixed Air 

Directional Airflow 

Displacement Ventilation 

Total supply aiiflow, L/s 

236 

236 

236 

Supply air particulate concentration, pg/ra^ 

0 

0 

0 

Supply air temperature, "C 

14 

14 

14 

Exhaust flow, Us 

236 

236 

236 

Exhaust temperature, "C 

31 

31 

31 

Particle concentration exhaust, pg/m^ 

140 

141 

140 

Particle concentration nonsmoking region, pg/m^ 

120 

98 

19 

Paidculate contaminant removal effectiveness (CRE) 

1.2 

1.4 

7.4 

CO concentration exhaust, pg/m^ (Ppm) 

1274 (1.1) 

1268 (1.1) 

1269 (1.1) 

CO concentration nonsmoking region, pg/m^ (Ppm) 

1081 (0.9) 

878 (O.S) 

157 (0.1) 

CO contaminant removal effectiveness (CRE) 

1.2 

1.4 

8.1 


Note: Bold indicates calculated results from input data. 
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(c) Displacement Ventilation 

Figure 2 Paniculate concentration fills at 1.7 m above fioor, 472 Us ventilation rate. 
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(a) Mixed Air 



(b) Directional Airflow 



(c) Displacement Ventilation 


Figure 3 CO concentration fills at 1.7 m above floor, 472 L/s ventilation rate. 
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Particulate concentration fills at 1.7 m above floor, 236 L/s ventilation rate. 
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The CRE values for directional airflow systems were 1.4 
at the lower ventilation rate and 2.7 at the higher ventilation 
rate. At the lower ventilation rate, the space.s had insufficient 
airflow to adequately control thediffusion of the contaminants 
from the smoking area into the nonsmoking area. Figures 2b 
and 3b display the concentration contours at the 472 L/s venti¬ 
lation rate, while Figure 4b shows tlie particulate coticeutra- 
tion contours for the lower supply air case. 

Displacement Ventilation 

The purpose of displacement ventilation is to introduce 
cool supply air into the space at low velocities, tlius minimiz¬ 
ing any air turbulence or mixing within the .space. I’he air w'as 
supplied close to the floor and tliereftire pooled along the floor. 
Air tlicn flowed vertically upward around all heat sources. 
ContaininantS from tobacco sources were heated in the 
combustion process and their natural tendency to rise as buoy¬ 
ant particles was enhanced by tlie vertical airflow pattern in 
displacement ventilation .spaces. 

The dispersion of contaminants into the non.smoking area 
was found to be minimal w'ith the use of displacement venti¬ 
lation. The highest CRE values (particulate and CO) resulted 
from displacement ventilation systems and were 7.4 at the 
lower ventilation rate and 17.5 at the higher ventilation rate 
(Table 1). Even operating at the lower ventilation rate, the 
displacement ventilation system provided nonsmoking area 
exposure levels that were lower than either the mixed air or 
directional airflow .system operating at their higher ventilation 
rate. Figures 2c, 3c, and 4c illusttate tlie typical concentration 
contours for this air distribution system at die 1.7 m elevation 
in the space. It can be observed that the spread of the particu¬ 
late and CO within the model was verj'similar at the two venti¬ 
lation rates. 

Movement of people in a space with displacement venti¬ 
lation may reduce the CRE value, since some air turbulence is 
created (Brohus and Kielsen 1994). It is not expected to 
change the relative differences in contaminant removal effec¬ 
tiveness for the three ventilation systems. A recent test room 
study where the effectiveness of mixing ventilation system 
arrangements was corapaied to displacement ventilation 
system arrangements using a smoking/nonsmoking zone set¬ 
up seems to conErm our findings (Enbom et al 2000). Enbom’s 
study found that some displacement system configurations 
would reduce nonsmokers’ exposure by 98% when compared 
with the situation of complete mixing. In our case, we calcu¬ 
late from Table 2 (4/63 xlOO) a reduction of 94%. 

CONCLUSIONS 

CFD modeling can be a useful tool to evaluate HVAC 
designs of spaces before they are built. It allows engineers to 
view air movement, room temperature, and concentration of 
contaminants within a space and to perform analyses of the 
impact of design changes. 
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The vast majority of all spaces employ the mixed air 
amangement. This study demonstrated that the other HVAC 
arrangements modeled, i.e., directional airtlow and displace¬ 
ment ventilation, may be alternatives for designing hospitality 
spaces where smoking occurs. In all cases modeled, the 
displacement ventilation system consistently had the highest 
CRE values (Table 1) and the directional airflow design 
consistently ranked above all mixed air systems. 

The diiectional airflow systems modeled showed CRE 
values more than twice as iiighas the mixed air systems at the 
higher ventilation rate. If a sufficient quantity of air is moving 
in one direction across the space, then contaminants generated 
in the smoking side of the space will be pushed toward the 
exhaust. The CRE of directional airflow systems is highly 
dependent upon the total airflow moving from the nonsmok¬ 
ing area into the smoking area. Higher airflow is expected to 
produce higher CREs, while at too low a flow, a result more 
like mixing should occur. Displacement ventilation performed 
best at directing contaminants to the upper regions of the space 
and was not as dependent upon total supply airflow as the 
mixed air and directional airflow systems were. 

Both particulate and gas phase movement were modeled 
with their different diffusion coefficients. There were small 
differences in the overall results between these substances as 
both contaminants moved in the direction of the prevalent air 
current. 

Usiug CRE is a straightforward method for determining 
and comparing how effective any ventilation system is in 
removing contaminants, or providing a partitioning effect 
within the space. This method should become a valuable lAQ 
measure in the future, as it i.s easy to calculate for any real 
space. This study demonstrated that ventilation rates affect the 
contaminant levels within a space, but siguificant improve¬ 
ments in minimizing concentration Of smoke in nonsmokiug 
areas can also be achieved by switching to either a diiectional 
airflow or displacement ventilation system. 
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